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The relationship between the nonperturbative Green's functions of Yang-Mills theory and the 
confinement potential is investigated. By rewriting the generating functional of quantum chro- 
modynamics in terms of a heavy quark mass expansion in Coulomb gauge, restricting to leading 
order in this expansion and considering only the two-point functions of the Yang-Mills sector, the 
rainbow-ladder approximation to the gap and Bethe-Salpeter equations is shown to be exact in this 
case and an analytic, nonperturbative solution is presented. It is found that there is a direct con- 
nection between the string tension and the temporal gluon propagator. Further, it is shown that for 
the 4-point quark correlation functions, only confined bound states of color-singlet quark-antiquark 
(meson) and quark-quark (baryon) pairs exist. 
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1. INTRODUCTION 

In quantum chromo dynamics [QCD] there are two primary aspects of the theory that remain to this day elusive: 
confinement and dynamical chiral symmetry breaking. Both involve the existence of a nontrivial, nonperturbative 
scale that is generated dynamically. Typically when one works with a physical theory such as QCD, one does not work 
directly with physical quantities, but rather with (the more readily calculable) gauge-fixed Green's functions. In the 
process of 'solving' the theory, one must regularize and rcnormalize at some scale, assigning the various parameters of 
the theory physical values such that the dynamically generated nonperturbative scale acquires physical meaning. This 
naturally raises the question: is there a simple relationship between the Green's functions and the nonperturbative 
physical scale determined externally? 

Further to this, Yang-Mills theory as the sub-theory of QCD without quarks does not have any directly observable 
quantities (the only observables are the glueballs and their spectrum, neither of which can be detected without having 
a detector made of matter, i.e., quarks and electrons). Yang-Mills theory is known to be confining in the sense that 
Wilson loops in lattice simulations exhibit an area law at large separations and this area law gives rise to a linearly 
rising potential, the coefficient of which is the string tension. Now whilst the lattice string tension exists, it has to 
be related to some physically observable scale in the real world and this has been explicitly done [1]. Unfortunately, 
the Wilson loop is not particularly amenable to calculation in continuum Yang-Mills theory (for recent work on this 
subject, see Ref. [2] and references therein), although it is the most natural quantity on the lattice. However, the 
physical interpretation of the Wilson loop of pure Yang-Mills theory as corresponding to two infinitely heavy quark 
sources connected via gluon exchange may provide a way of sidestepping the difficulties of considering the (gauge 
invariant) Wilson loop in the continuum and providing the connection between the (gauge-fixed) Green's functions 
of Yang-Mills theory and physical quark confinement. The heavy quark mass expansion (and the effective theory 
associated with it) is a well established area of QCD and indeed uses phenomenologically motivated potentials in 
place of the Yang-Mills sector [3-5] (see also, for example, Ref. [6]). Moreover, in the heavy quark mass limit the spin 
of the quark decouples from the system (the so-called heavy quark symmetry) leading to dramatic simplifications. 

In full QCD, the Bethe-Salpeter equation is one tool for describing systems composed of quark pairs with interactions 
given via the Green's functions of the theory. In the light quark sector where one is concerned primarily with dynamical 
chiral symmetry breaking and the spectrum of the light mesons, covariant gauge Bethe-Salpeter studies have proven 
to be extremely useful for phenomenology (e.g., Refs. [7, 8] or the contemporary review Ref. [9] and references therein). 
Such studies are typically restricted to the rainbow-ladder truncation, where the Bethe-Salpeter kernel is reduced to 
the single exchange of a dressed gluon and so that the problem is tractable. However, it is known that the Bethe- 
Salpeter kernel is not so simple and indeed, there has been recently some considerable attention focused on going 
beyond the rainbow-ladder approximation both in terms of vertex corrections (see for example Refs. [10-15]) and in 
unquenching effects (e.g., Refs. [11, 16-18]). Whilst more sophisticated kernels do entail a considerable increase in 
effort, one of the more intriguing results of these studies is that (aside from the meson decay mechanism induced by 
unquenching [17]) the rainbow-ladder truncation seems surprisingly robust. In Coulomb gauge, there also exist studies 
of dynamical chiral symmetry breaking and the Bethe-Salpeter equation (see for example Ref. [19] or the more recent 
Ref. [20] and references therein) . Here, potential models (corresponding to the rainbow-ladder truncation) have been 
used. 

As was recognized early on. Coulomb gauge is an ideal choice for studying nonperturbative phenomena [21], but 
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because of the inherent noncovariance, significant technical progress has been only recently achieved (mainly in Yang- 
Mills theory). The Gribov-Zwanziger confinement scenario [22-24] is particularly relevant in Coulomb gauge. In 
this scenario, the temporal component of the gluon propagator is infrared enhanced (providing a long-range confining 
force) whereas the transverse spatial components are infrared suppressed (and therefore do not appear as asymptotic 
states). The relevance of Coulomb gauge stems from the fact that the system reduces naturally to physical degrees 
of freedom [24]. Using functional methods, it has been formally shown that there exists a conserved and vanishing 
total color charge (along with the absence of the infamous Coulomb gauge energy divergences) [25]. In addition, 
the Dyson-Schwinger equations have been explicitly derived [26-28], along with the Slavnov-Taylor identities (for 
Yang-Mills theory) [29] and one-loop perturbativc results have been calculated [27, 28, 30]. On the lattice, there 
are also initial results for the Yang-Mills propagators now available (although at present the lattice sizes are still 
modest) [31, 32]. Of note is that on the lattice, the temporal gluon propagator appears largely independent of the 
energy (in the noncovariant Coulomb gauge, the propagators are in general dependent on both the energy and spatial 
momentum as separate variables, unlike in covariant gauges) and is consistent with a l/iT^-behavior in the infrared. 
Furhermore in Coulomb gauge the static spatial lattice gluon propagator is found to vanish in the infrared and to 
agree with Gribov's formula [31]. The lattice results support the findings of the variational approach to continuum 
Yang-Mills theory in Coulomb gauge [33, 34] (c.f. also [35, 36]). 

To establish a connection between the Green's functions of Yang— Mills theory and the physical world represented 
by the existence of the string tension that confines quarks, we propose to study here the full QCD Bethe-Salpeter 
equation with a heavy quark mass expansion at leading order in Coulomb gauge. Given that all quarks are confined 
irrespective of their mass, it is reasonable to suppose that the existence of the confinement potential and the value 
of the string tension are independent of the mass and the configuration of the quarks, the only important quantity 
being their separation in this respect. 

Now, the heavy quark mass expansion is clearly not a new concept, nor is the Bethc-Salpcter equation and the 
existence of a confinement potential. We shall be working in Coulomb gauge and again, there do exist studies of 
the Bcthe Salpetcr equation within this context as discussed above. Thus, we should be specific about what is new 
in this work. Heavy quark effective theory is done typically in covariant gauges. Separately, in Coulomb gauge the 
study of the Bethe-Salpeter equation has focused on chiral symmetry breaking and the existence of the confinement 
potential has been taken as being more or less a scittlcd matter (see for example Ref. [19]). However, the Coulomb 
gauge Bethe-Salpeter studies have not gone beyond the leading (rainbow-ladder) approximation. In this study, we 
shall work with heavy quarks in Coulomb gauge and with the emergent results for the explicit Green's functions of 
the Yang-Mills sector (at least under truncation). In addition we will be working nonpcrturbativcly, deriving and 
utilizing in part the Slavnov-Taylor identity for the quark-gluon vertex and its peculiar Coulomb gauge form. We thus 
combine the separate areas of study. What we shall show is that as a result of this combination, the leading (rainbow- 
ladder) approximation is exact in the case considered, leading directly to the interpretation of quarks being confined 
by a linearly rising potential and providing an explicit link between the Green's functions of the nonpcrturbative 
Yang-Mills sector and the external physical scale. 

The paper is organized as follows. In Section 2 we consider the generating functional of Coulomb gauge QCD at 
leading order in the heavy quark mass expansion. The tree-level quark and antiquark propagators are derived and 
the relevant Feynman prescriptions are introduced. In Section 3 the nonpcrturbative quark gap equation is discussed. 
Using the Slavnov-Taylor identities in Coulomb gauge (derived explicitly in Appendix A) and truncating the Yang- 
Mills sector to include only dressed two-point functions, it is found that the rainbow approximation to the quark 
and antiquark gap equation is nonpcrturbativcly exact. This result is confirmed by the semipcrturbativc approach 
considered in Section 4. In Sections 5 & 6, the Bethe-Salpeter equations for quark-antiquark and diquark states are 
considered. It is shown that the ladder approximation to the Bethe-Salpeter kernel is nonperturbatively exact and 
that only color-singlet meson and SU{2) baryon states have finite energy solutions. Moreover, a direct connection 
between the temporal gluon propagator and the string tension is established, at least under the truncation scheme 
considered. 



2. HEAVY QUARK MASS EXPANSION 

In this study, we are concerned with combining nonperturbative physics with the heavy quark mass expansion. The 
two are not automatically compatible; for example, whilst the mass expansion assumes that the mass parameter is 
the largest scale in the problem, loop integrals involve momenta up to the UV-cutoff scale (supposing for the moment 
that a UV-cutoff regularization is employed) which is also assumed to the the largest scale. In heavy quark effective 
theory [HQET] , the recognition of this apparent contradiction leads to the so-called matching procedure [3] . This is 
just one example of the effect of the heavy mass expansion and some care must be taken in proceeding. The overriding 
concern here is that we wish to use wherever possible the full nonperturbative QCD functional formalism, i.e., the 
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gap and Bethe-Salpeter equations and Slavnov- Taylor identities for the complete quark fields, rather than HQET 
expressions that refer to the heavy quark degrees of freedom. The effect will be (see below) that we are obliged to 
restrict ourselves to leading order in the mass expansion. This naturally prohibits a quantitative analysis involving 
real quarks (which are not infinitely heavy) but does not deflect from the goal of studying the relationship between 
the string tension and the Yang-Mills Green's functions. 

Throughout this work, we shall use the conventions and notations established previously in Refs. [26-29]. We work 
in Minkowski space with the metric gfj,^, = diag(l, —1) and all minus signs associated with the spatial components of 
covariant and contravariant vectors are explicitly extracted. Roman sub/superscripts (z, j,. . . ) denote spatial indices, 
superscripts (a, b,.. .) denote color indices in the adjoint representation. Configuration space coordinates may be 
denoted with subscript {x, y,. ■ ■) when no confusion arises. Other index notation will be explained as the context 
arises. The Dirac 7-matrices satisfy {7^,7"} = 2g^''. The explicit quark contribution to the full QCD generating 
functional is [28] 

Z[X,X] = J I'^expjj J d^xq^ix) ij°Dq + ij-D - m ^^qpix) 

X exp |z J d'^x \Xoc{x)qa{x) +q^{x)xa{x)] + iSym^ (2.1) 
with the temporal and spatial components of the covariant derivative (in the fundamental color representation) given 

by 

Do = do-tgT"c7%x), 

D = V + z5TM"(a;). (2.2) 

In the above, generically denotes the functional integration measure over all fields present. qa{x) denotes the full 

quark field {q is the conjugate or antiquark field) where in this case, the subscript a refers to the fundamental color, 
fiavor and spin attributes collectively. The sources Xaix) and Xai^) are for the complete quark fields. The Yang-Mills 
contribution to the generating functional (which will only be of direct relevance to this study in Appendix A) is in the 
standard, second order formalism [27, 29]. Additionally, A and a refer to the spatial and temporal components of the 
gluon field, respectively. The T" are the (Hermitian) generators of the SU{Nc) group satisfying [T",T''] = if^-^'^T'' 
(with the fully antisymmetric structure constants /"'"') and normalized via Tr(T''T^) = S""^ /2. For later use we 
introduce the color factor associated with the quark self-energy: 

N'^ - 1 

ft = ^. (2.3) 

Now consider the following decomposition of the quark and antiquark fields: 

fc(x) = e-*™"" [h{x) + H{x)]^ , K{x) = e'"-o [P+g(x)]„ , H^{x) = e™^« [P-g(x)]„ 
q^{x) = e^--" \h{x)+Hix)]^, K{x) = e-""-o [q{x)P+]^ , H„{x) = e"^™-" [q{x)P-], (2.4) 

where the (spinor) projection operators are 

P± = i(l±70), P++P_=l, P+P_=0, Pl=P±. (2.5) 

This decomposition is a particular case of the heavy quark transform underlying HQET [3] . There, one recognizes that 
a heavy quark within a hadron is almost on-shell and moves with the hadron velocity v such that the 4-momentum 
can be written p^^ = mw^ + where |fc| ^ m|?;| and = 1 (such that when \k\ = 0, = m^). One then uses the 
general projectors P± = (l±)(^)/2 and with the exponential terms e^"""'^. The case used here corresponds to the rest 
frame of the quark, w** = (1,0), but within the context of the generating functional is simply a choice of (arbitrary) 
decomposition that will prove useful in Coulomb gauge. In fact, this choice will result in the simplification whereby 
the spatial components of the Yang-Mills Green's functions are absent at leading order in the mass expansion. The 
virtue of the heavy quark decomposition is that the projection operators satisfy the following further relations 

P+7°P+=P+P+, P+7°P-=0, P+yP+=0 (2.6) 

such that the following relations hold for the components of the quark field: 

hj°h = hh, H'y°H = -HH, hj°H = H-f°h = hfh = HYH = 0. (2.7) 
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Inserting the decomposition of the quark fields given by Eq. (2.4) into the generating functional Eq. (2.1) and using 
these relationships one obtains 

Z%x] = J I>$exp|i J (fx K{x)[iDo]^^h^{x) 

+K{x)\i')-D\ Hi3{x) +Hoc{x)\i^-D\ hi3{x)+Ha{x)[-2m-iDo]^gHij{x) 

L J a/3 L J a/3 

: I d^x [e-""-°x„(a;) [h{x) + H{x)]^ + e^^^'^ \h{x) +H{x)\^Xo.{x)\ +iSym^- (2.8) 



X exp < I 



The difference between the mass expansion (as used here) and HQET can now be explained. Because of our insistence 
on retaining the source terms for the full quark fields, our generating functional has only been rewritten in terms of 
different integration variables (note that the Jacobian of the transform is field independent and thus trivial) and not 
altered, although the source term expression in the above is slightly modified by the appearance of the exponential 
factors. We will thus be able to work with the full gap and Bethe-Salpeter equations for quarks which are derived 
from Z[x,x]- This is in contrast to HQET where the quark sources are replaced with sources for the /i-ficlds (the 
if-fields arc integrated out as below). The differences are mainly cosmetic at the level of this study (which will 
concentrate on the leading order in the mass expansion, see below) but serve to illustrate some useful points which 
will be highlighted as they arise. 

The point of the heavy quark decomposition of the fields is that for the /i-field ('large') components, the quark mass 
parameter m does not appear directly. When we integrate out the if-fields however, we get the following expression 

z[x,x] = 



V^Bet [iDq + 2m] exp |z J d'^x ha{x) [iDq]^^ h/3{x) 
h{x)i^-D + e-'™^°x(a;)l [iDq + 2m] J \i^-Dh{x) + e""^°x(a;) 

; j d^X [e-'"'''°Xaix)ha{x) + e'"'-''h„{x)Xa{x)] +tSYMy 



+ 



X exp < I 



(2.9) 



Obviously, since we have integrated out a nontrivial component of the original quark field, our expression is nonlocal 
and this is where the heavy mass expansion is necessary. We notice that [3] 



Det[«£>o + 2m] ~ O {l/m'^) , 
[iDo + 2m]-'^X0 ~ ±.x^ + 0{l/m'). 



(2.10) 



Leaving the e^'^"^^° factors as they are, we can thus write 



Z[x,x] = J V^expL J d^x 



K{x) [iDo]^. hffix) + — h{x)i^-D + e-'^^"x{x) il-Dh{x) + e'™^«x(a;) 



exp |i J d^x [e-'"'^'>Xc,{x)hc,{x) + e™^°/i„(x)xa(x)] + iSym^ + O (l/m^) . 



(2.11) 



Our generating functional is now local in the fields and arranged in an expansion in the parameter 1/m (this will 
be referred to as the mass expansion although strictly speaking, it is an expansion in the inverse mass). However, 
locality does not mean that the above expression can be directly applied. Let us consider the classical (full) quark 
field in the presence of sources: 



j V^qa{x) exp{iS} 



1 SZ 

Z 6iXa{x) 



e-""^''/i„(a;) + 



2m 



i^-Dh{x) + e'"'''°x{x) \ exp{iS} + O (l/m^) . 



(2.12) 



One sees immediately that even at 0(l/m), the classical quark field has components that involve interaction type 
terms {Dh ~ Ah) brought about by the truncation of the nonlocality. Of course, this is nothing more than the 
statement that the /i-field is nontrivially (and dynamically) related to the full (/-field. It also means that if we want 
to use the nonperturbative gap and Bethe-Salpeter equations (i.e., those equations derived from the action for the 
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full quark fields q and their sources x ^-nd which are the equations of QCD as opposed to HQET) then wc cannot 
expect that the mass expansion can realistically be extended far beyond the leading order in order to do practical 
calculations. As stated previously though, the aim is to investigate the connection between the Yang Mills sector and 
the physical world made of quarks; the string tension that represents our goal is not dependent on the quark mass 
(both light and heavy quarks are confined in the same way, as far as we know). Therefore, we restrict our attention 
to the leading order in the mass expansion as follows (and writing Do explicitly): 

Z[X,X] = Jv^expl^iJ d''xha{x)[ido,+gT''a''{x)]^0hi^{x)^ 

X exp |^ J d^x [e-'"'^°Xaix)hc{x) + e'"^°7^„(a;)x«(a;)] + iSym^ + O (1/m) . (2.13) 

The standard machinery of functional methods is now employed. From the observation that up to boundary terms 
(which are assumed to vanish) the integral of a derivative vanishes, we have the quark field equation of motion: 

0=1 exp{«<S} 

J Slha {x) 

= I V^[[ido, + 5TV«(a;)]„^ hp{x) + e'^^OxaW} exp{i<S} + O (1/m) . (2.14) 



The field equation of motion for the antiquark gives equivalent results and can be neglected. The generating functional 
of connected Green's functions (for the full quarks) is defined via Z = e^. Denoting the derivatives with respect to 
sources of W[x, x] with an obvious bracket notation, the classical fields are: 

qc{x) = 1?$9„(x)exp05} = =:<iXa{x)> 
Also {p is the source for the cr- field, implicit within Z and W), 
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V<!>a^{x)q^{x) exp {iS} = . = Z[x, x] [<^/o"(x)^x„(x)> + <ip'^ {x)><iXa{x)>] • (2.16) 

dip''(x)dlXa{x) 



Notice that from Eq. (2.13), the following relations hold: 

SZ[x,x] _ f^^^-rmxoi. r„ci , rntM^\ ^Z[x,X. 



I D$e-»™^o/i„(x) exp {iS} + O (1/m) , = " / '^'^e""''°ha{x) exp {iS} + O (1/m) (2.17) 



(recalling the earlier discussion of neglecting the O (1/m) terms) such that the field equation of motion, Eq. (2.14), 
can be written in terms of derivatives of W: 

= b5o.]„^e^™^° <iX0{x)>+[gTXpe''"'"' [<ip''{x)iX0ix)> + <ip'' {x)><iXis{x)>] + e^'"^<'Xa(^) 

+e»(l/m). (2.18) 

Factoring out the exponential terms gives then 

= [idox - m]^p <iX0{x)> + [<ip''{x)iX0{x)> + Kip'' {x)><iX0{x)>] + Xa{x) + O (1/m) . 

(2.19) 

To continue, wc make a Lcgendre transform in order to construct the effective action for the full quark fields. Explicitly 
separating the quark and Yang-Mills components: 

r[^,q,q] = W[J,x,x] - ^Ja^a - iXaQc, - iq^Xoc- (2.20) 

such that 

1oi{x) =<lXa{x)>, Xa{x) = - <iq^{x)>, 
Qa{x) = - <lXa{x)>, Xa{x) =<iqa{x)>, 

G^ix) =<ip''{x)>, p"'{x) = - <ia''{x)> (2.21) 
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where the same bracket notation for derivatives of the effective action with respect to the classical fields is used (there 
is no confusion between the two sets of brackets since the two sets of derivatives are never mixed) . Jq, and denote 
generic gluonic (Yang-Mills) sources and classical fields, respectively; the index a here refers to all attributes of the 
(gluonic) object in question (including its type and position argument) and we use the common convention that this 
index is either summed or integrated over as appropriate. The field equation of motion can then be rewritten in terms 
of derivatives of the effective action (which are the proper Green's functions when sources are set to zero): 

<iq^{x)>= [ido^ - m]^^ qp{x) + [gT"]^^ [<ip''{x)iXp{x)> +a%x)qp{x)\ + O (1/m) . (2.22) 

With the field equation of motion written in the above forms, we can now derive the Feynman rules for the quark 

components of the theory (the Yang-Mills parts are already known [27]). We start with the quark propagator. From 
Eq. (2.19), ignoring interaction terms and functionally differentiating we get the tree-level propagator in configuration 
space 

= - m]^f, <^X7(^)*X^(2;)>^°' -i5^J{z - x) + O (1/m) . (2.23) 

We now write the Fourier transforms to define both the anti- and quark propagators, respectively, utilizing the 
translational invariance: 

<tX0{^)tXa{x)>=Wgg0^{z - x) = J dke-'''<'--^Wgg0^{k) (2.24) 
such that in momentum space, we get 

= J dke-^''<^-^^{[-ko-m]^^W^^^^{k)-z6,^} + 0{l/m), 

= J cike-^'<^-y^ {w^''X{k)[ko-m]^^+iSf^,} +0{l/m). (2.25) 
Naively, one would write the solutions as 

However, in order to define the Fourier transform, one must first define a prescription for handling the poles in the 
energy integral. Unlike the conventional tree-level quark propagator, we do not have a pair of simple poles in the 
complex fco-plane, we have instead single poles and this is due to the mass expansion. For the quark propagator, we 
write 

Notice the following. The tree-level quark propagator is a scalar quantity (or rather, diagonal in the outer product 

of the fundamental color, flavor and spinor spaces) as a consequence of the mass expansion. This is one of the most 
striking features of HQET and physically corresponds to the decoupling of the spin from the heavy quark system. 
W^'^^ is identical to the heavy quark tree-level propagator [3] up to the appearance of the mass term. In HQET, where 
one uses sources for the /i-fields directly, this term docs not appear and is simply due to a shift of the energy by an 
amount m (such that one is working with the 'small' momenta around the rest mass energy). We retain this term 
for completeness in order to see what effect it may have - the Yang-Mills sector may not be compatible with such a 
shift in the energy, although this will turn out not to be the case. Also, note that the kinetic term of the tree-level 
propagator would read — fc^/2m in the denominator factor and is at higher order in the mass expansion. Such terms 
are obviously important to the UV properties of the loop integrals but will not play any role here since we shall be 
interested in the infrared limit. 

Returning to the issue of the Feynman prescription, one of the direct consequences of the single energy pole is 
that when one considers a closed quark loop (a virtual quark- ant iquark pair) at lowest order, the energy integral 
automatically vanishes (see also Ref. [3] for an alternate discussion on this topic), i.e., 

" 0. (2.28) 



[ko — m + le] [ko +po — m + le 
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It can be justified that all such closed quark loops vanish in similar fashion (indeed, precisely this type of integral will 
be heavily used in the following Sections); in other words, the heavy mass expansion at leading order is quenched. 
One interpretation of the above result is that due to the presence of only the single pole, the heavy quark can only 
propagate forward in time and the corresponding antiquark is not present. This leads to a conundrum: how can 
we construct a quark-antiquark pair in the Bethe-Salpeter equation? The resolution is straightforward - in the 
Bethe-Salpeter equation one is not considering a virtual quark-antiquark pair but rather a system composed of two 
separate unphysical particles. It is known that in Coulomb gauge. Gauss' law forbids the creation of a colored state 
in isolation (the total color charge is conserved and vanishing [25]) and so, the Feynman prescription for the quark 
(or the antiquark) propagator has no physical meaning in isolation. The quenching of the theory means that closed 
quark loops vanish and from covariant gauge Bethe-Salpeter studies, the connection between such loops and the 
Bethe-Salpeter kernel is qualitatively understood [17]. In the quenched case, the quark and antiquark lines of the 
Bethe-Salpeter equation are never connected by a primitive vertex (unlike the closed quark loop). Another way of 
seeing this is to consider the flavor nonsinglet Bethe-Salpeter equation - in other words the quark and antiquark 
constituents considered as two distinct flavors. Because the quark and antiquark arc genuinely separate (unphysical) 
particles, we are at liberty to choose that the antiquark propagator in the Bethe-Salpeter equation has the opposite 
Feynman prescription as follows: 

As shall be seen, this results in a physical interpretation for the quark-antiquark Bethe-Salpeter equation as a whole. 

Another way of viewing the above is to realise that the four components of the full quark Dirac equation (quark 
and antiquark moving either forwards or backwards in time according to causality) have been separated by the 
mass expansion (the 'large' h or 'small' H field components) and the breaking of the time reversal symmetry. In the 
context of the mass expansion, there is only the quark or antiquark moving forward in time corresponding to the above 
Feynman prescriptions. The closed quark loop involves a quark going backwards in time (similarly for an antiquark 
loop) which is prohibited, or more precisely suppressed by the mass, so that such closed quark loop integrals vanish 
at leading order whereas the quark-antiquark Bethe-Salpeter equation has a physical solution. 

For the proper two-point and three-point functions, we use functional derivatives of Eq. (2.22) and we get directly 
in momentum space 

^fleik) =t[ko- m] + O (1/m) , rg^„^(A:) = i[ko+ m] S^p + O (1/m) , 
rSl%(^i' ^3) = [gTX^ + O (1/m) , rg);„^(fci, h) = - [gT^]^^ + O (1/m) . (2.30) 
Importantly, the tree-level spatial quark-gluon vertex does not appear at leading order in the mass expansion: 

4X/3i-r(Jl,, = 0(l/m). (2.31) 

Notice the ordering of the indices for the Fg^^ vertex. Also notice that ri°|^^(fc) = ~r^^^^^(— A;) since the two-point 
function requires no Feynman prescription and is diagonal in the outer product of the nmdamental color, flavor and 
spinor spaces. In addition, WqqTqq = 1 as usual. 

Now, because of our insistence of using the full quark sources, all the nonperturbative equations involving the 
quarks (the Dyson-Schwinger equations for two-point and three-point functions, the Slavnov-Taylor identities and 
the Bethe-Salpeter equation) will not alter their form at leading order in the mass expansion (beyond leading order, 
as previously illustrated, the nontrivial classical fleld and the resolution of the Legendre transform would play a role). 
The only alterations are that the tree-level factors have been changed and we demand that only the leading order 
terms in the mass expansion of the resulting equations are retained. Obviously, in order to solve the nonperturbative 
system we must further specify our truncation scheme. Recall that at leading order in the mass expansion, the theory 
is already quenched. The truncation scheme we propose is to consider only the dressed two-point functions of the 
Yang-Mills sector (i.e., the nonperturbative gluon propagators derived from a hypothetical solution of the complete 
Yang-Mills theory) . This amounts to setting all the pure Yang-Mills vertices and higher n-point functions occurring 
in the quark equations to zero. It is worth pointing out that because of the mass expansion, any loop diagram with 
a tree-level spatial quark-gluon vertex will be suppressed, such that the number of loop diagrams arising because of 
the Yang-Mills vertices is heavily restricted and those loops that do contribute will not include the leading order 
perturbative corrections (the fully temporal gluon Green's functions Ta-a-aT^aa-aa, ■ ■ ■ are zero at tree- level). The most 
important physical point that will emerge is that when we set the Yang-Mills vertices to zero, we exclude the non- 
Abelian part of the charge screening mechanism of the quark color charge and any potential glueball states. On the 
other hand, the charge screening mechanism and glueball contributions of the gluon field (i.e., the color string) is 
implicitly encoded in the nonperturbative form of the temporal gluon propagator. 
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With the truncation scheme as outhncd above, the Yang-Mills sector collapses to the inclusion of a single object: 
the temporal gluon propagator which is written as [27] 

W:',{k) = S'^'^D^^k'). (2.32) 

There are three important features to this propagator. Firstly, there are indications on the lattice that the dressing 

function D„a is largely independent of energy [32], justifying the energy independence of the above form. Indeed 
on general grounds, the temporal gluon propagator must have some part that is constant in the energy in order to 
cancel closed ghost loops and resolve the Coulomb gauge energy divergences (see also the more formal considerations 
of Ref. [37]). Second, the lattice analysis indicates that the dressing hmction Dg-o- is infrared divergent and is likely 
to behave as for vanishing k^. We are interested mainly in the relationship between _D(j(j (as the inpiit of the 
Yang-Mills sector) and the string tension so we will not need the specific form until towards the end. Third, the 
product g^Daa is a renormalization group invariant quantity and thus a natural candidate for being relevant to the 
physical string tension [24, 29]. 



3. GAP EQUATION: NONPERTURBATIVE TREATMENT 

Let us begin by considering the Dyson-Schwinger equation for the quark two-point proper function (the gap 
equation). In full QCD (i.e., second order formalism, Coulomb gauge without the mass expansion and derived from 
the first order formalism results of Ref. [28]), it reads {duj = rf'*w/(27r)^): 

T-,,c.6{k) = r£)„,(fc) {r(°);„^(fc, - fc)H%,^^(a;)r|,,^,(a;, -fc, k - u:)W:l{k - 

+rfi%^(^' ^ - k)W^qP,{^)^'^gA-,Sj{^, -k, k - ^)W%\,^{k - ^)] (3.1) 

(Waa is the spatial gluon propagator, which will be unimportant here). The full quark-gluon vertices obey the 
Slavnov-Taylor identity. The derivation of this identity follows from the invariance of the action under a Gauss-BRST 

transform [24] that is peculiar to Coulomb gauge and is similar to the identities derived for the pure Yang-Mills sector 
[29]. This derivation is presented in Appendix A (and is a technical result of this study in its own right). The identity 
reads, Eq. (A.18): 

km,,^p{k,MM) = l^TI^Aafi,{k„k2,ks)T^^{-k^) 

ks 



+Tqgas{kl) f|;7;cq(fcl + 90, - 9o; ^2) + igT'^ 

+ 



S0 



T%^cqik2 + qo,k3-qo;ki)-igT^ rggSf3i-k2) (3.2) 

where ki + k2 + k^ = 0, qo is the arbitrary energy injection scale, F-^ is the ghost proper two-point function, F^-ccg 
and Tq-ccq are ghost-quark kernels associated with the Gauss-BRST transform (see Appendix for details). 

In order to use the Slavnov-Taylor identity as input for solving the gap equation, we first apply our truncation 
scheme in the context of the heavy mass expansion at leading order. Starting with the dressed spatial quark-gluon 
vertex, consider the terms that contribute to the Dyson-Schwinger equation shown schematically in Fig. 1. According 
to the truncation scheme, we set all Yang-Mills vertices to zero, meaning that diagrams (c-f) are excluded. This 
then leaves us with the tree- level term (a) and the quark- loop term (b). However, both of these involve at least one 
tree-level spatial quark-gluon vertex which is not present at leading order in the mass expansion. Thus, we obtain 
the nonperturbative result that 

T^qAa0i{ki,k2,ks) = O{l/m). (3.3) 

It is straightforward to justify that the ghost-quark kernels of the Slavnov-Taylor identity, given their definition, 
Eq. (A. 19), involve Yang Mills vertices and do not contribute here. Thus, in our truncation scheme and at leading 
order in the mass expansion, the Slavnov-Taylor identity reads 

k°srlg,„0{ki,k2, fcs) = Tqq^siki) [igT"] - [igT"] F5g,^(-fc2) + O (1/m) . (3.4) 

Clearly, the truncation scheme thus results in an Abelian type Ward identity. Moreover, since the temporal quark- 
gluon vertex is simply multiplied by the temporal gluon energy (the primary feature of Coulomb, as opposed to 
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(a) (b) (c) (d) 




(e) (f) 



FIG. 1: Diagrams that contribute to the Dyson-Schwingcr equation for the spatial quark-gluon vertex (without prefactors 
or signs). Internal propagators are fully dressed and blobs represent dressed proper vertex and (reducible) kernels. Internal 
propagators represented by springs may be either spatial (A) or temporal {a) propagators, deished lines represent the ghost 
propagator and solid lines represent the quark propagator. 

covariant, gauge Slavnov-Taylor identities) and the quark proper two-point function is color diagonal, we can imme- 
diately write the solution: 

^,aa0{h, k2,h) = ^ {T-^ [T^,{h) - r^,(-fc2)]}„^ + O (1/m) . (3.5) 

The above solution is trivially satisfied at tree-level. We notice however, that there appears to be a potential problem 
with the energy. When ^3 = 0, but fcs ^ (a spacelike gluon configuration), the simple pole must somehow be canceled 
by the difference of proper quark two-point functions. Since the spatial momentum configuration is arbitrary, this 
means that Tqq{k) — >■ Tgq{ko) + (l/m). The demand that the nonperturbative vertex solution to the Coulomb gauge 
Slavnov-Taylor identity be free of kinematic divergences (here, simply the l/fcg factor) is a variation of the familiar 
covariant gauge situation considered in Rcf. [38]. 

Inserting the results, Eq. (3.5) and Eq. (3.3), for the vertices, using the Feynman rules given by Eq. (2.30), with the 
temporal gluon propagator given by Eq. (2.32) and resolving the color structure the nonperturbative gap equation, 
Eq. (3.1) under truncation and at leading order in the mass expansion thus reads 

^gqasiko) =l[ko- m] SaS - Q^Cp [ ^'^^""^^ 't^Wqqapi'^o) \Pqq{'^o) " ^qq{ko)\Bi + O (1/m) . (3.6) 

J {ko - uJo){k - ujy 

There exists one particularly simple solution to this equation as we shall now demonstrate. It is given by 

Wqqaf){k)= ~' +0(l/m), rqq„p{k)=i6c.0[ko-C] + O{l/m) (3.7) 

and where the constant {duj = d?uj/{2'K)^) 

C = m+ \g'C, I ^^^^ + O (1/m) . (3.8) 

Putting the above solution into the Slavnov-Taylor identity, we also have that for the vertex 

rlq,a.^{kuk2, fcs) = [gT'] ^^ + (1/m) . (3.9) 
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In other words, the dressed temporal quark-gluon vertex is trivial and the gap equation reduces to the rainbow 
truncation. There is a subtlety to the solution, Eq. (3.7), involving the ordering of the limits in the spatial and 
temporal integrals and potential divergences in the constant given by Eq. (3.8). Here, let us consider the case when 
we perform the temporal integral first under the condition that the spatial integral is somehow regularized and finite 
(this will be done throughout the rest of the paper). Inserting the solution given by Eq. (3.7), the gap equation under 
truncation, Eq. (3.6), and with a regularized (i.e., finite) spatial integral denoted by the subscript r, reads 

r , 2r< f ^f^Daaik-u) i dujo , ^nri/ \ 

Cr = m + q Cf / =; — ^ t:— hm / 7 r+C 1/m) 

^^^^ .CF[^^^4^ + Oil/m). (3.10) 

^ Jr {k — cjy 

The effect of performing the temporal integral first is that the regularized constant Cr in the denominator factor 
becomes irrelevant. The spatial integral now involves no external scale and removing the spatial regularization, we 
arrive at our above result, Eq. (3.8). 

A brief discussion about the physical interpretation of these results is in order. Firstly, it might be the case that 
there exist other solutions to the truncated gap equation. As will be shown in the next Section though, the above 
solution can also be derived from a semi-perturbative type of expansion. However, as is the case with many systems, 
especially those with strong couplings or phase transitions, the fully nonperturbative solution might not be the same 
as the resummed perturbative solution. Second, that the solution involves potentially divergent constants is not a 
comfortable situation but does not necessarily contradict the physics. We notice that the quark propagator has a 
single pole, so cannot represent physical propagation (which requires a covariant doiiblc pole) and this arises obvioiisly 
from the truncation of the mass expansion and where the charge conjugation symmetry has been explicitly broken. 
Indeed, the position of the pole has no physical meaning since the quark can never be on-shcU. That the single pole 
is shifted to infinity simply means that either one requires infinite energy to create a quark from the vacuum or that 
should one have an incoming quark (from some other hadron), only the relative energy is important. Note that in 
Coulomb gauge it is known that the total color charge of the system is conserved and vanishing [25] , so one cannot 
prepare an isolated colored state of a single quark under any circumstances. Further, the divergences here have no 
interpretation with regards to renormalization, at least within the context of the mass expansion to leading order. 
The mass parameter cannot be renormalized simply becaiisc one cannot construct an appropriate counterterm in the 
action that is linear in m to absorb the constant C. Also, the quark field renormalization is trivial at leading order, 
as one sees from the explicit form of the temporal quark-gluon vertex, Eq. (3.9). The upshot of this is that one must 
consider only the relative energies in the system the divergence of the absolute energy has no physical meaning and 
this will become especially relevant when we study the Bethe Salpetcr equation (indeed this has been known for quite 
some time, see for example Ref. [19]). 

Having discussed the quark propagator, let us now discuss the antiquark propagator. Recall that at tree-level, we 
used a different Feynman prescription for the two denominator factors and this gives rise to some rather interesting 
physical consequences. As previously discussed, the heavy mass expansion employed here breaks the charge con- 
jugation symmetry relating particle and antiparticle, so we cannot expect that the two propagators are necessarily 
equivalent. Starting with the gap equation for full QCD, Eq. (3.1), we reverse the ordering of the quark and antiquark 
functional derivatives that form the quark Green's functions (still within the context of the full quark fields and 
sources) and rearrange the ordering to get the gap equation for the antiquark propagator: 

-r,^Sa{-k) = -rgL(-^) - / ^'^ {i^%.Sj{-k, OJ, k - a;)W,^^^(-a;)rg);^ J-a;, k,oj - k)Wf,{k - w) 

+^%A5yj{-k,uJ, k - uj)W,-,^^{-uj)T^^20aii-^^k,uj- k)Wt,,^{k -co)]. 

(3.11) 

Applying our truncation scheme reduces the above to 

%5a(-fc) = rg)j-fc) + 1 di^T''^^,s-y{-k,co,k-uj)W,^^0{-u;)r'^g^p^^^ (3.12) 
In similar fashion, we have the Slavnov-Taylor identity for the antiquark-gluon vertex: 

-klr%.0o.ik2,k„k3) = +r,^0s{k2) [igT'^j'^sa - [wT%r,^Sa{-h) + 0{l/m) . (3.13) 
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(a) (b) 



FIG. 2: Diagrams that contribute under our truncation scheme to the Dyson Schwinger equation for the temporal quark-gluon 
vertex (without prefactors or signs). Internal propagators are fully dressed and blobs represent dressed proper vertex and 
(reducible) kernels. Internal propagators represented by solid lines represent the quark propagator. 

The form of the solution to Eq. (3.12) is similar to the previous results: 

W,-,^p{k) = + O (l/m) , r,5«/3(^) = i5ap [ko+Cr] + O (1/m) (3.14) 

[fco + Cr + ie\ 

with the corresponding solution for the vertex 

T%^^^{kiM, ks) = - [gT'']^^ + O (l/m) . (3.15) 

However, the constant (obtained in the same way as before by first doing the temporal integral under the assumption 
that the spatial integral is regularized in some hypothetical manner) is now 

Cr = m- \g-C^ [ i^^^ + O (l/m) . (3.16) 

^ Jr ^ 

We notice that the sign of the loop correction has changed and this will turn out (in the context of the Bethe-Salpeter 

equations for mesons and diquarks) to have far-reaching consequences. For the moment wc interpret this rcisult as 
being simply another manifestation of the breaking of the charge conjugation (particle-antiparticle) symmetry via 
the mass expansion. 



4. GAP EQUATION: SEMIPERTURBATIVE TREATMENT 

Whilst we have already solved the gap (and anti-gap) equation, it proves instructive to reconsider it within the 
context of a semiperturbative analysis since this will naturally introduce a technical feature crucial for considering 
the Bethe-Salpeter equation nonperturbatively. We have seen that the solutions for the proper two-point function 
leads to a temporal quark-gluon (and antiquark-gluon) vertex that is not dressed. Under our truncation scheme, 
the nonperturbative Dyson-Schwinger equation for the temporal quark-gluon vertex involves the diagrams shown in 
Fig. 2. The semiperturbative expansion is based on a hybrid loop expansion whereby all internal propagators are 
taken to be dressed, but all internal vertices are tree-level. To show that all the loop corrections (contained in diagram 
(b) of Fig. 2) vanish, it suffices to consider two types of diagram, given in Figs. 3 and 4. 

In Fig. 3 we consider a single ladder exchange correction to the temporal quark-gluon vertex. This diagram 
(neglecting the overall color and prefactors) gives rise to the following scalar integral (as in the previous section, the 
spatial integral is regularized such that one can perform the temporal integral first): 

f S uj D„„{k - J) 1 /"^ (kcii 

Jr {k - a;)2 2n J_^ [wq - Cr + te] [wq + qo-Cr + ts]' 

Now we apply the following identity (for finite, real a, b; the case a = 6 is trivial): 



r°° dz 1 r°° f 1 1 \ 

J-oo [z - a + is][z -b + ts] {a - b) \[z-a + ie] [z-b + ts]j 



0. (4.2) 
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FIG. 3: Ladder type loop correction to the temporal quark-gluon vertex. Internal propagators are fully dressed: solid lines 
represent the quark propagator and springs denote the temporal gluon propagator. 




FIG. 4: Generic crossed box type loop correction to the temporal quaxk-gluon vertex. Internal propagators are fully dressed: 
solid lines represent the quark propagator and springs denote the temporal gluon propagator. The box represents any combi- 
nation of interactions allowed under our truncation. 



Thus we see that where the spatial integral is regularized, the temporal integral vanishes. It is simple to see that 
the planar one-loop diagrams with two or more external temporal gluon legs (which under the truncation scheme 
considered here connect only to the internal quark line) and one internal temporal gluon will also vanish. 

Now let us consider a generic crossed box (nonplanar) type of diagram, illustrated in Fig. 4. Considering only the 
temporal double integral components of the explicit internal quark propagators, we have the following form 



J — c 



[vo - ai + te] [wo -a2 + le] [ojq -h go - as + [^o -vo-po-a4 + le] 

duja 

^ [ujo ~ a2 + «e] [wo + qo - as + le] [uq - po - ai - + 2ie] 



I dvo \ 



-27rz 



ai + le] [vo - Wo + Po + 04 - . 

dtuo 



[cjo -a2 + ts] [u>o -h go - as + «e] [wo - po - ai - 04 -h 2ie] 







(4.3) 



where in the last line, we have used a variation of the identity Eq. (4.2). Thus we have the result that the generic 
crossed box typo of diagram shown in Fig. 4 also vanishes. 

Given that both the single ladder type exchange diagram and the generic crossed box diagrams considered so far 
vanish, it is easy to see that any vertex dressing diagram will vanish (including all subdiagrams such as internal vertex 
corrections and so on), since all diagrams are merely variations or combinations of these two under our truncation 
scheme. This result is a consequence of the fact that the energy and Feynman prescription of the denominator factors 
follow the quark line through the diagram so that eventually the identity, Eq. (4.2), can be used. It is also precisely 
the reason why all closed quark loops vanish, as previously discussed. 
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FIG. 5: Homogeneous Bethe-Salpeter equation for quark-antiquark bound states. Internal propagators are fully dressed and 

solid lines represent the quark propagator. The box represents the Bethe Salpeter kernel K and filled blobs represent the 
Bethe-Salpeter vertex function V with the (external) bound state leg given by a dashed line. See text for details. 

Thus, the semiperturbative expansion confirms our previous result that the temporal quark-gluon vertex remains 

bare to all orders. Also elear should be that the result applies to the antiquark-gluon vertex too. With the corre- 
sponding simple forms for the self-energy integrals of the gap and anti-gap equations (which as we recall, reduce to 
the rainbow truncation), the results for the quark propagator functions are also confirmed. Notice though that whilst 
the all orders semiperturbative result must match the nonperturbative result, the converse is not necessarily true. It 
remains the case that there may exist further solutions, but these must be purely nonperturbative in character if they 
exist. That the identities that have been introduced here conform to the nonperturbative case previously studied is 
useful since it allows us to apply them with confidence to the Bethe-Salpeter equation in the next Sections. 



As has been emphasized, because we initially consider the generating functional for full quark sources and fields 
we are at liberty to use the full functional (nonperturbative) equations as a starting point and then subsequently 
apply our mass expansion and truncation scheme. Let us now consider the homogeneous Bethe-Salpeter equation for 
quark-antiquark bound states. In full QCD and with our conventions (the minus sign arises from the definitions of 
the Legendre transform and Green's functions), this equation reads 



where fc+ = fc + ^P, k- = + (^ — l)f (similarly for p±) are the momenta of the quarks, P is the pole 4-momentum of 

the bound state (assuming that a solution exists), ^ = [0, 1] is the so-called momentum sharing fraction that dictates 
how much of the total meson momentum is carried by each quark constituent, K represents the Bethe-Salpeter kernel 
and r is the Bethe-Salpeter vertex function for the particular bound state that one is considering and whose indices 
explicitly denote only its quark content (we will see later what color, flavor and spin structure the solutions may 
have). Physically, results should be independent of ^ and this has been numerically observed in phenomenological 
studies [7]. The Bethe Salpeter equation is shown pictorially in Fig. 5. 

Aside from the quark propagators (which we shall discuss shortly), the central element to solving the Bethe-Salpeter 
equation involves the construction of the kernel K. As discussed in the Introduction, for technical reasons the most 
widely studied system is based on the ladder kernel which is either constructed via the interchange of a single gluon 
(for example [8]) or as a phenomenological potential (see for example Ref. [19]). However, there has been much 
recent attention focused on the construction of more sophisticated kernels. One key element of the construction is 
the axialvector Ward-Takahashi identity [AXWTI], which relates the gap equation to the Bethe-Salpeter kernel and 
which ensures that chiral symmetry and its spontaneous breaking are consistently implemented (e.g., Refs. [10, 15, 19]). 
Here, we shall show that the ladder Bethe-Salpeter kernel is exact at leading order in the heavy mass expansion and 
under our truncation scheme. This derivation follows in the same way as the semiperturbative analysis of the previous 
Section. 

Following from the truncation of the heavy mass expansion to leading order, it was arguc;cl that the antiquark 
propagator must be treated as distinct from the quark propagator. This means that we have to be very explicit about 
which propagator is which. Since we are studying the quark-antiquark system (we will analyze the quark-quark, or 
diquark system in the next Section), the Bethe-Salpeter equation for our purposes more properly reads 



5. BETHE-SALPETER EQUATION: MESONS 




dkK^p,s-y{p, k; P) [Wgg{k+)T{k; P)Wgg{k-)l 



(5.1) 




dkK^p.Sjip, k; P) [W^g{k+)T{k; P){-l)W^^{-k-)] 



(5.2) 
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FIG. 6: Generic crossed box type of diagram that contributes to the Bethe-Salpeter kernel. Internal propagators are fully 
dressed, whereas vertices are tree-level. The upper (lower) solid line denotes the quark (antiquark) propagator; springs denote 
the temporal gluon propagator and the box represents any combination of nontrivial interactions allowed under our truncation 
scheme. See text for details. 



where we have explicitly identified the antiqiiark propagator contribution (it corresponds to the lower line of Fig. 5) by 
reordering the functional derivatives and there is implicitly a similar antiquark contribution within the kernel which 
absorbs the explicit minus sign. Recall that we are implicitly considering only the flavor non-singlet case. 

To construct the full Bethe-Salpeter kernel, consider the generic semiperturbative crossed box contribution given 
in Fig. 6. As before, at leading order in the mass expansion and under our truncation scheme we only have the 
temporal quark-gluon and antiquark-gluon vertices, both of which have been shown to be given by their tree-level 
forms. Additionally, it was justified that all planar diagrams with multiple external temporal gluon legs vanish. Thus, 
aside from the ladder contribution to the Bethe-Salpeter kernel, this generic crossed box diagram contains all possible 
nontrivial contributions. Such a diagram has at least the following terms in the; temporal intcigral (as before, we 
assume that the spatial integral is regularized and finite so that we are able to firstly perform the temporal integral 
without complication) 



/ 



(5.3) 



[uo +p%-Cr + is\ ... [ojo -k^+Cr+ te] ■ 



The first factor corresponds to the explicit quark (upper) propagator, the last factor to the explicit antiquark (lower) 
propagator. Implicitly (represented by the dots), there may be multiple propagator factors which carry the same 
dependence on the integration energy wq but crucially, all these will have the same relative sign for the Feynman 
prescription term, i.e., wq + is, regardless of whether they originate from internal quark or antiquark propagators. 
Therefore, this type of integral can always be rcdiiccd to the difference of integrals over a simple pole and with the 
same sign for carrying out the analytic integration, just as in Eq. (4.2) and the semiperturbative treatment of the 
vertex discussed in the last Section. Thus, all generic crossed box diagrams in the Bethe-Salpeter kernel are zero 
and one is left with simply the ladder contribution to the kernel. Actually, this could have been anticipated from 
the beginning - the AXWTI connects the self-energy term of the gap equation and the Bethe-Salpeter kernel and 
since it has been shown explicitly that the self-energy integral reduces to the rainbow truncation, the corresponding 
Bethe-Salpeter kernel is simply given by ladder exchange. 

Given all this, the Bethe-Salpeter equation for the quark-antiquark system, at leading order in the mass expansion 
and within our truncation scheme, can be explicitly written as 

+0{l/m). (5.4) 

Inserting the nonpcrturbative results for the propagators and vertices so far, Eqs. (3.7,3.14,3.9, 3.15) and taking the 
form, Eq. (2.32), for the temporal gluon propagator, we get the equation (the temporal and spatial integrals are 
separated as before) 

We sec immediately that the flavor and spin structure of the meson decouples from the problem - this is well a known 
property of the heavy mass expansion. The color structure will be discussed shortly. Since the external energy po 
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docs not enter the right-hand side of the above equation, we further have that the Bethe-Salpetcr equation must be 
independent of the relative quark energy (and only implicitly dependent on the putative bound state energy Pq) . 
Thus we can write 



ra/3(p;-P) = g 



2 / dkD^^ip- k) 



{p- 

ip-k) 



TT{k;P)T'' 



dko 



Jr (V- k' 



ai3 2n J_^ [k^ -Cr + is] - - is] 
T'T(fc; P)T'' 



+ 0{l/m) 



010 



{p-kf [Po-Cr+Cr + 2ls] 



0(l/m). 



(5.6) 



Thus, at leading order in the mass expansion, inserting the expressions Eq. (3.10) and Eq. (3.16), for Cr and Cr, 
respectively, it is now clear that 



diU Daa{L0) 



r„0(p;P) = -<?2 



dkDaa{p- k) 

{p- ky 



TT{k;P)T'' +0{l/m). 



(5.7) 



We notice that the explicit quark mass contributions of the self-energy expressions cancel. This is a feature of the 
quark-antiquark Bethe-Salpeter equation - it does not make any reference to the origins of its constituents and why 
for example the pion can be a massless bound state of massive constituents. Physically, one can visualize that the 
quark and antiquark are moving with equal and opposite velocities such that the centre of mass system (the bound 
state) is stationary. This is related explicitly to the choice of Feynman prescription for the constituent quark and 
antiquark. Were the Feynman prescription for the antiquark chosen to coincide with that of the quark, the right-hand 
side of Eq. (5.6) would simply vanish and there would be certainly no physical quark-antiquark state. The Feynman 
prescription for the antiquark corresponds precisely to a particle moving with the opposite velocity. Also, at leading 
order in the mass expansion, the momentum sharing parameter, ^, has dropped out. Thus, the results retain the 
important physical requirement that they be independent of ^. Shifting integration momenta, we can write 

PoT^M P)=g^ I '^'^^r^'^^ {CfTo^M P) - [T^np- P)T«]„ J + O (1/m) . (5.8) 

J r ^ 

To see the physical meaning of this equation, we rewrite the Bethe-Salpeter vertex function as a Fourier transform: 

r„0{p;P)= J dye-'P-yrc0{y) (5.9) 

(in the homogeneous Bethe-Salpeter equation, the total momentum P denotes the solution and is not a variable). 
We also write for the color structure 



[T''mTXf3 = CMr.^(y) 

where Cm is yet to be identified. Then, the Bethe-Salpeter equation reduces to 



(5.10) 



J dye-'P yPoTc,^{ 
with the simple solution 



= J dye-'P-yy 



di0Daa{uj) 



{CfTo^M - e^-yCMT^p{y)] + O (1/m) (5.11) 



^wL'crcr(w) 



{CF-e'^-yCM} + 0{llm). 



(5.12) 



Because the total color charge of the system is conserved and vanishing [25], neither the quark or antiquark can 
exist as an independent asymptotic physical state. Thus, the bound state energy, Pq, can only increase linearly as 
the separation between the them increases (physically confining) or be infinite when the hypothetical regularization 
is removed (so that the system cannot be physically created). Whether the system is confining or disallowed can 
only depend on the color structure, since the temporal gluon propagator dressing function would be common to both 
situations. In configuration space, an infrared confining solution is characterized by the solution Pq = o'\y\ for large \y\ 
and where a is the string tension such that as the separation between the quark and antiquark increases, the energy 
of the system should increase linearly without bound and infinite energy input is required to fully separate them (at 
least in the absence of unquenching). The small |^ (and large |a;|) properties are of no concern here. The Fourier 
transform integral needed for the infrared confining solution is 



J Stt' 



(5.13) 



16 



If the temporal gluon propagator dressing function is more infrared divergent than l/|a5|, as above, then 

Cf = Cm (5.14) 

is required such that the spatial integral is convergent and the energy of the system well-defined, as the hypothetical 
regularization is removed (since we are interested in the low |aJ| regime, it becomes clear that the regularization here 
would be infrared in character). Using the Fierz identity for the generators T": 

2 = ^«7'5^/3 - ^^aph^ (5.15) 



gives the condition 



CpTa-iiy) = CM^ajiV) = ^SajTpisiy) - ^Tajiv), (5.16) 



or with the definition Eq. (2.3), 

^aiiy) = Sa'jl'iy). (5.17) 

In other words, the quark-antiquark Bethe-Salpeter equation with D^a more infrared singular than l/\oj\ can only 
have a finite solution for color singlet states where the divergent constant integral coming from the unphysical quark 
self-energy cancels; otherwise the energy of the system is divergent. 

Another way to see that only color-singlet states are physical is to consider that if the temporal gluon propagator 
is the origin of a potential in configuration space, this potential can always be shifted by some spatial constant [19]. 
In momentum space, this means that one can make the replacement 

g2Daa^ ^ g^£^^ + const X {27rfS{oj) (5.18) 

without changing the physical bound state energy. This automatically gives the constraint Cp = Cm regardless of 
whether or not the temporal gluon propagator is infrared enhanced. 

Assuming that in the infrared (as is indicated by the lattice data [32] or by the above argument about the non- 
existence of asymptotic quark states), D^^ = XjuP' where X is some combination of constants (and further knowing 
that g^X is a renormalization group invariant [24, 29]), then 

Po = c7|y| = ^^^|y|+0(l/m). (5.19) 

The above result is that there exists a direct connection between the string tension and the nonperturbative Yang- 
Mills sector of QCD at least under the truncation scheme considered here. The veracity of the truncation scheme will 
be discussed at the end. 



6. BETHE-SALPETER EQUATION: DIQUARKS 

Let us now briefly take a look at the diquark Bethe-Salpeter equation. The difference between this and the 
previously considered quark-antiquark system is actually rather simple from a technical standpoint, but leads to a 
completely different physical result. Recall that the quark and the antiquark propagators share the same Feynman 
prescription relative to their energy and only the constant components change. Thus, preserving the momentum 
routing of the Bethe-Salpeter equation, the result that the crossed box contributions to the Bethe-Salpeter kernel 
extends to the diquark case since this is purely dependent on the Feynman prescription. This means that we can 
immediately write down the Bethe-Salpeter equation for diquarks, at leading order in the mass expansion and within 
our truncation scheme: 

+0{\lm). (6.1) 

Again, the indices of the Bethe-Salpeter vertex function correspond to the quark content of the diquark and since 
the flavor and spin content decouple from the system, we shall only be interested in the color content of the diquark 
below. Expanding this out as before, we get the analogous result 

Po - 2m - g-Cp [ ^"^-^"^ 1 T„M P) = g' f [T^^, [T%^ V^P- ^l^ + O (1/m) . (6.2) 

J r ^ J J r ^ 
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Fourier transforming as previously, and writing 

[TXx r,«(^) = CnT^M (6-3) 

gives the solution 

Po = 2m + g^ [ ^'^^^^'^^ {Cf + e^^ yCu) + O (1/m) . (6.4) 

The reappearance of the quark mass simply indicates that in contrast to the quark-antiquark system, there are now 

two co-moving quarks. The expression for the anti-diquark system is identical to the above, but with minus twice 
the mass - their velocities are simply reversed. The diquark is antisymmetric under interchange of the two quark legs 
and since the spin and flavor structure decouples, this means that the color structure must be antisymmetric. As for 
the quark-antiquark system, the system can only have a finite energy confining (energy increasing with separation) 
solution (if Cd = —Cp) or no finite solution at all. Using the Fierz identity Eq. (5.15), the color condition then reads 

-CfTc^M = CdTc^M = \vp^{y) - i^TaM- (6-5) 

Demanding the diquark color antisymmetry and with the deflnition Eq. (2.3) this becomes 

Nj - Nc-2 = Nc = -1, 2. (6.6) 

So in SU{Nc = 2) there exists a confined, antisymmetric bound state of two quarks - the SU{2) baryon and otherwise 
there are no (finite) physical states. 

7. SUMMARY, DISCUSSION AND CONCLUSIONS 

The connection between the Green's functions of Coulomb gauge Yang-Mills theory and physical quark confinement 
has been studied. Since all quarks are confined (irrespective of their mass), the generating functional of full QCD 
was expanded in the mass parameter and the leading order considered in a manner appropriate to Coulomb gauge 
and such that the system simplifies dramatically. This allowed the usage of the full nonperturbative quark equations 
of QCD in order to study the confinement properties. Using the Coulomb gauge Slavnov-Taylor identity for the 
quark-gluon vertices (whose derivation is one of the results of this study) and truncating the Yang -Mills sector to 
include only the nonperturbative gluon propagator, it was shown how the rainbow approximation to the quark and 
antiquark gap equations is exact in this case. It was then demonstrated that the corresponding ladder approximation 
to the Bethe-Salpeter equation was also exact. The mass expansion breaks the time-reversal and charge conjugation 
properties of the theory and in order to describe physical quark-antiquark states, the relevant Feynman prescription 
was introduced. With the analytic solutions to the gap equation, the; Bethe-Salpeter equation was solved for the 
quark-antiquark and quark-quark channels. It was found that the only solutions correspond to confinement, namely 
that only color-singlet meson and SU{2) baryon states have finite energy. This energy must increase linearly with 
separation cither because: (i) there are no asymptotic quark states as is known in Coulomb gauge from the conserved 
and vanishing total charge [25] , or (ii) as is suggested from the lattice results [32] , the temporal propagator is infrared 
enhanced. Further, there exists a direct connection between the temporal gluon propagator of Yang-Mills theory and 
the string tension, at least within this truncation scheme. 

Because of the (almost embarrassingly) simple results, it is worth re-emphasizing the input to the calculations. 
Different areas of study have been combined here: the heavy quark mass expansion, the Dyson-Schwingcr formalism 
along with the Slavnov-Taylor identity and the Bethe-Salpeter equation, all of which have been considered within 
Coulomb gauge. Indeed, the simplification arising from this combination underscores how powerful the choice of 
Coulomb gauge is for nonperturbative studies. There are two approximations used. Firstly, the temporal gluon 
propagator is considered to be energy-independent. This is consistent with available lattice data [32], although 
explicitly in disagreement with perturbation theory [27]. Since this study is interested in the infrared properties 
of the theory, using the lattice suggestion seems justified. Actually, the connection between the explicit one-loop 
perturbative expressions and the leading order perturbative expansion (in covariant gauges, the resummation of 
the leading logarithms via the anomalous dimensions) is not yet understood in Coulomb gauge. Further, since the 
nonperturbative ghost propagator in Coulomb gauge is strictly independent of the energy [29] in order to cancel the 
closed ghost-loops occurring in the Yang-Mills expressions, the temporal gluon propagator must have some part that 
is also independent of the energy (the spatial gluon propagator is explicitly dependent on the energy). Thus it seems 
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reasonable to assume that whatever energy dependence the temporal gluon propagator does have (and that is not 
visible with current lattice configurations), it is not important to the conclusions here. 

The second approximation used in this study was to neglect the vertex and higher n-point functions of Yang Mills 
theory. Recall that the tree-level spatial quark- gluon vertex is suppressed by the mass. This means that the number 
of truncated terms at leading order in the mass expansion is heavily suppressed. It was seen that without such terms, 
the quark-gluon vertices (and the Bethe-Salpeter kernel) reduced to their tree-level forms. This suggests that whilst 
the dressing function Dcro- does implicitly contain all nonperturbative effects associated with the dynamical dressing 
of the color charge (including, for example, potential glueball states), the quark-gluon vertices correspond to the 
naked quark color charge. One might express the results here as corresponding to a dressed color string confining 
two naked color sources. In the gap and Bethe-Salpeter equations, the explicit Yang-Mills contributions presumably 
cannot obliterate the effect of the rainbow-ladder exchange for a given temporal propagator input since (assuming 
the infrared enhancement as suggested from the lattice [32]) one would require an exact 'anti-confining' cancellation. 
Covariant gauge studies of the effect of non-Abelian corrections to the Bethe Salpeter equation as applied to the light 
quark sector indicate that they are dominant compared to the Abelian corrections [11]. However, in the case studied 
here, there are no Abelian corrections and the diagrammatic content of the non-Abelian contribution is reduced by the 
mass expansion, making a quantitative comparison to known Bethe-Salpeter studies somewhat speculative. Perhaps 
a better comparison about the nature of neglecting the non-Abelian corrections is with Ref. [24], which refers to the 
Wilson loop of pure Yang-Mills theory. With effectively the same truncation (i.e., neglecting the Yang-Mills vertices), 
it was demonstrated that a string tension (the so-called Coulomb string tension) can also be extracted. Subsequently, 
it was shown that the physical string tension is lower than the Coulomb string tension [39] (a statement colloquially 
referred to as "no confinement without Coulomb confinement"). Thus, we anticipate that the effect of including the 
non-Abelian corrections to the formalism presented here would result not in the removal of the linearly rising bound 
state energy (i.e., not the cancellation of the ladder exchange) but rather in the softening of the coefficient by shifting 
the pole position by some finite amount. Physically, this amounts to the charge screening of the quark color charge. 

The outlook of this study is extremely positive, with several avenues that might prove fruitful for further inves- 
tigation. Maintaining the truncation scheme, one can consider the three quark (or two quark plus antiquark etc.) 
Faddeev equation to investigate SU (3) baryons and whether or not the simple results for color confinement persist to 
the many-body system. Even more speculatively, given the simplicity of the leading order mass expansion one might 
imagine going even further and considering the two quark, two antiquark system. A second important line of study 
would be to consider the explicit inclusion of the Yang-Mills vertices to the combined system of Dyson-Schwinger, 
Slavnov-Taylor and Bethe Salpeter equations, cither phenomenologically or self-consistently to assess the influence 
of the truncation scheme considered so far. A third possibility is to investigate what happens when one chooses 
a different velocity parameter in the original decomposition of the quark fields. This would mix the temporal and 
spatial parts of the Coulomb gauge framework and allow one to investigate the connection between the respective 
gluon propagators. 
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The derivation of the Slavnov-Taylor identity for the quark-gluon vertex will be presented in this Appendix. This 

derivation follows almost automatically as an extension to the Slavnov Taylor identities for the Yang-Mills sector 
[29]. The full QCD action in the standard, second order formalism reads 
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Appendix A: Quark Slavnov-Taylor identity 




(A.1) 



where the (antisymmetric) field strength tensor F is given in terms of the gauge field A^: 



(A.2) 




(A.3) 
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such that for infinitesimal 9%, the fields transform as (recall that a = Aq 

„ 1 



ahc—hnc 



9a 



-doe" - 

9 

9; = 9a-t^°[n„^g^, 



Fixing to Coulomb gauge (V-A° = 0) via the Faddeev-Popov technique introduces new terms into the action: 



where 



D 



ab ca6 



(A.4) 



(A.5) 



(A.6) 



is the spatial covariant derivative (in the adjoint representation), A" is a Lagrange multpher field to locally implement 
the gauge condition, c" and are the Grassmann-valued ghost fields. The action is invariant under a Gauss-BRST 
transform [24] whereby the infinitesimal spacetime-depcndcnt parameter 0^ is factorized into two Grassmann-valued 
components: 0'^ = cJJJAj, where SXt is the time- dependent infinitesimal variation. The Gauss-BRST transform is a 
Coulomb gauge specific form of the standard BRS transform, allowed because the gauge-fixing term does not involve 
temporal differential operators. The variations of the new fields read: 



1 



r'"'c''c'=SXt, SX" = 0. 



Including a source term. 



d^x 



5, = / 

the generating functional is given by 



Z\J]= j D$ exp {iSqcd + iSfp +iSs}. 



(A.7) 



(A.8) 



(A.9) 



The Coulomb gauge Slavnov-Taylor identities arise from regarding the Gauss-BRST transform as a change of inte- 
gration variables under which the generating functional is invariant. Since the Jacobian factor is trivial [26] and only 
the source term varies, we deduce that 



-I 



5 [i5Xt 



exp {iSqcd + iSpp + iSs + i5Ss} 



j V^e^{iSQCD + tSpp + iSs) j d^xSit - X,) |-i {Oy,) c« + /"^^^c^ 

TsaX 



g-^ix^ixCx ~^ f '^ix^ix'^x ^Xax^x'^a, 



x-ipx^ paA.ax "T X^TJ^ o ^x'-'x^x 

9 ^ 



(A.IO) 



Just as for the pure Yang-Mills case, the time dependence of the variation, SXt , results in the (5-function constraint 
S{t — Xo)- The above identity is most usefully expressed in terms of proper functions and repeating the manipulations 
of Ref. [29] , one arrives at the identity 



= / d^xS{t -xo) I i (as <<>) ct - r"' <<> [<^P^^rx> +<^lcl] - ^ 



(-VS) 



<iAt> 



-r'"= <iA1^> t,,{x) [<iJ%vfn> +A)^cl] - -XI <c> <icl> [<^^^^r?> +^0^] 

9 ^ 
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Note that functional derivatives involving the Lagrange multiplier result merely in a trivial identity such that the 
classical field can be set to zero [29]. Further, to derive the quark Slavnov-Taylor identities, one functional 
derivative with respect to id^ is needed and then the ghost fields/sources can be set to zero. Implementing this then 
gives 



= / d^x6{t -xo) \-- {d° <iai>) 5{z - x) - /«*^ <i<> 



<ipli%> -ialb'"'5{z - x) 



1 

+- 

9 



<iAt^> 



-lT^0 <iqax> 



^ <^WTfx> +5''''5{z - x)iq0. 



I3a 



<iXlSxiV> -S'^'Siz - x)iqp 



<nax> 



} 



(A.12) 



Two further functional derivatives with respect to iq.y^ and iqg^ are taken and all remaing fields/sources set to zero. 
Given that [29] 



= 0, 



(A.13) 



j=o 



one obtains the Slavnov-Taylor identity for the quark-gluon vertices in configuration space: 



=Jd^xS{t-xo) 



<tqs^iq^uiai>) S{z - x) + 



1 



<i^icf> 



+iT^p <iq5yiqocx> 

(52 



> 



5iqsJici 



<iXi3xrnx> 



+ S'^'Siz - x)S-y^d{io - x) 



J=0 



6'^"-diz - x)6sfid{v - x) 



(A.14) 



As with all Slavnov-Taylor identities, the above involves nontrivial ghost scattering-like kernels. Introducing some 
notation, these can be expressed via the Legendre transform in terms of a loop integral over proper Green's functions 
and propagators in the following way: 



y{x,z,uj) = igT^p 



Siq-yu^id^ 



<^X/3x^Vx> 



j=o 



WTap {- <^X/3x^Xe><^qe^Q'1Ul^'^'x><Jx^JK><ni^VT><^CT^cil^K> 
+ <^X|3x'^X^i><^lK'^<ll^'^CT^ci><lrf^lr]r>} , 
5^ 



diqgjicf 



<^X0xtVx> 



{<r]lir]r><lCrtc'^t^K><lJKtJx><iqSvn£i^X><iXe^Xl3x> 

- <ir]lir]r><iCri4iqsyiq^><ix^iX0x>} igT^^- 



(A.15) 



In the above, J and $ are dummy sources/fields refering to either A or cr only, the internal indices otherwise refer to 
all attributes of the object in question (summed or integrated over). The configuration space Slavnov-Taylor identity 
can thus be written 



= / d'^xS{t - xq) 



<iqsyiq-yoji(Ti>) S{z - x) 

<iq5v'<-<lax> 



<iqsyiq^ujtA'}^> 



+ 



^rccqa-iix, z, w) + igT^^S{z - x)S{u} - x) 

^-.ccqSai^^ 2' ^) - WTsJ{z - x)5{v - x) <«g„x«57U'>} • 

Defining the Fourier transform for the vertex functions (all momenta incoming): 

T{x, y,z) = j ah dk2 dk3 (27r)^^(fci +k2 + fc3)e-''=i-^-'*=^-^-''=^-^r(A;i, fcs, fcs) 



(A.16) 



(A.17) 
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one can write the Slavnov-Taylor identity in momentum space (dropping the dirac and fundamental color indices for 

convenience) 

+r^,(fci) [f^;-,,(fci + go, h - go, k2) + igT'''] + [f + go, ^3 - go, h) ~ igT'''] r^,(-fc2). (A.18) 

In the above, ki + k2 + = 0, Fee is the proper ghost two-point function and go is the (arbitrary) energy injection 
scale that arises from the time-dependence of the Gauss-BRST transform. One can see the strong similarities of this 
expression with the Yang-Mills case [29]. In momentum space, the kernels can be written 

X [r|j?g,(a;,p2,Pi - w,p3) -r|^«(w,P2,-P2 - w)W^«a(P2 +w)r|,;,(pi -W,P3,P2 +w)] , 

f g;ccq(Pi>P2,P3) = j duj [r^^^(a;,p2, -P2 - uj)W^l{p2 + u;)T^g^{p3,pi - a;,p2 + w) - T^igg{uj,p2,P3,Pi - w)] 

xWl\oj)Wgg{iJ - Pl)^flT^ (A.19) 

where the indices k, A, refer to the gluonic field types tr or A (with the associated spatial index). Just as in the 
Yang Mills case, the above Slavnov Taylor identity, Eq. (A.18), in conjunction with the kernels can be solved in 
principle to give the temporal quark-gluon vertex Tgg^ in terms of purely spatial, ghost or quark propagators and 
proper functions. The identity is trivially satisfied at tree-level (using the Feynman rules of Ref. [28]) and can be 
verified at one-loop perturbatively (it is a purely technical exercise, so not suitable for presentation here). 
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